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Abstract 

Q\ . The Crab Nebula has been observed by the HEGRA stereoscopic system of 4 imaging 
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air Cerenkov telescopes (IACTs) for about 100 hours from 1997 September to 1998 
March. The recent detailed studies on the system performance give an energy thresh- 
old and an energy resolution for 7-rays of 500 GeV and ~ 18%, respectively. The 
Crab spectrum was measured with the IACTs system in a very broad energy range 
from 500 GeV to 25 TeV, using the observations at zenith angles up to 50 degree. 
The pure power-law spectrum in TeV 7-rays constrains the physics parameters of the 
nebula environment as well as the models of photon emission. 

1. Introduction 

The Crab Nebula has been observed and studied over an exclusively broad photon 
energy range embracing radio, optical, X-ray bands as well as high energy 7-rays up to 
hundreds of TeVs. The various theoretical scenarios of photon emission are primarily 
based on the Synchrotron Compton model, which combines the synchrotron and 
inverse Compton emissions from high-energy electrons, accelerated up to ~ 100 TeV, 
which interact with magnetic field and seed photons within the nebula (for review 
see Aharonian & Atoyan, 1998). The predicted spectrum in the TeV energy domain 
appears to be very sensitive to the model parameters: the value of magnetic field, the 
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nature of seed photons, the maximum energy of electrons, etc. Thus a measurement 
of TeV Crab spectrum sets major constrains on theoretical expectations. 
The imaging air Cerenkov technique was used for observations in the TeV energy 
range (for a review see Weekes et al, 1997). Currently the Crab serves as a standard 
candle for TeV 7-ray emission. However the Crab energy spectrum, as measured 
by several groups, still differs in power-law slope as well as in absolute normaliza- 
tion. While using the same fundamental principles, presently operating imaging air 
Cerenkov telescopes have a different energy threshold (e.g. CAT: ~ 200 GeV; Whip- 
ple: - 250 GeV; HEGRA: - 500 GeV; CANGAROO (large zenith angles):- 7 TeV) 
and dynamic energy range. In addition the energy spectrum measurements can be 
affected by the systematic errors, limiting the accuracy of the result. Stereo imaging 
gives several advantages for the spectrum studies, comparing with the single telescope: 
(a) direct measurement of shower impact parameter; (b) good energy resolution; (c) 
wide dynamic range; (d) extended abilities for systematic studies using several im- 
ages for an individual shower. The HEGRA IACTs system was primarily designed 
for detailed spectrum measurements in the TeV energy domain utilizing the advan- 
tages of the stereoscopic observations. Note that detailed systematic studies for the 
spectrum evaluation technique have been recently made using Mrk 501 observations 
in 1997. Here we present the Crab data for the last observation period, analysed with 
the recently developed technique of energy spectrum evaluation. 
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2. Observations 

The Crab observations reported here were made with the stereoscopic system of 4 
HEGRA imaging air Cerenkov telescopes, which are located at La Palma, Canary 
Islands (Daum et al., 1997). Each of the 4 telescopes consists of a 8.5 m 2 reflec- 
tor focussing onto a photomultiplier tube camera. The number of photomultipliers 
in the camera was 271, which were arranged in a hexagonal matrix covering a field 
of view with a radius of 2.°3. Any telescope camera was triggered when the sig- 
nal in two nextneigbour of the 271 photomultiplier tubes exceeded a threshold of 8 
photoelectrons, and the system readout started when at least two telescopes were 
triggered by Cerenkov light from air shower. The detection rate was 12.6 Hz near the 
zenith. The Crab was observed in a wobble mode; i.e., the telescopes were pointed in 



Table I. Observation time 



N 


Period 




Number of tel-s 


Number of runs 


Time, hr 


1 


1997 Sep - 


1998 Jan 


4 


152 


45.63 


2 


1997 Oct - 


Nov 


3 


48 


14.26 


3 


1998 Feb - 


March 


4 


84 


22.62 



Declination ±0.5° aside from the nominal Crab position (a sign of angular shift was 
altered from one run of 20 min to another). This is useful for continuous monitoring 
of the cosmic-ray background taking the OFF-source region being symmetric about 
the camera center, and 1° apart relative to the ON-source region. Observations of 
the Crab at zenith angles up to 50 degree were made from 1997 September 1 to 1998 
March 29, for a total of 82.5 hr of data taken at good weather. By the fire accident 
at the HEGRA site, one of the telescopes was damaged and was out of operation 
for a month in October-November 1997. At that time, Crab observations were made 
with only three telescopes in the system, providing an event rate of 10.3 Hz near the 
zenith. Due to unstable weather and substantial amount of dust which came from 
Sahara desert to the island, the average detection rate in 1998 February-March was 
reduced down to 10.7 Hz in observations near the zenith. The total exposure times 
for the three periods are summarized in Table 1. 

3. Analysis and result 

The stereoscopic imaging analysis of the data is based on the geometrical reconstruc- 
tion of the shower arrival direction and the shower core position in the observation 
plane, as well as on the joint parametrization of the shape of the Cerenkov light im- 
ages. The simultaneous registration of several (>2) Cerenkov light images from air 
shower provided an angular resolution of ~ 0.°1 for 7-ray showers. The first 7-ray se- 
lection criterium used here is 2 < 0.1 [deg 2 ], where O 2 is the squared angular distance 
of the reconstructed shower arrival direction to the source position. This orientational 
cut is very loose and gives a 90% 7-ray acceptance whereas the number of cosmic ray 
showers is reduced by a factor of about 55. In addition we analysed the data by the 
"mean scaled width" parameter, < w >. This parameter was introduced in order to 
provide an almost constant 7-ray acceptance over the dynamic energy range of the 
telescope system. Thus, the second 7-ray selection criterium was < w > < 1.4, which 
accepts most of the 7-ray showers (~ 99%) while the relative acceptance of cosmic ray 
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showers is 26%. The calculation of parameter < w > includes the zenith angle depen- 
dence of the image shape. The collection function of energy and zenith an- 



Table II. Crab data summary. 



Period: 


Z.A. 


Time, hr 


N 7 


N cr 


R 7 , hr" 1 


R cr ,hr 1 


S/N* 




up to 25° 


24.9 


2691 


5493 


107 


220 


23 


1 


25° - 40° 


14.2 


929 


3276 


65 


230 


11 




40° - 50° 


6.5 


360 


1355 


55 


208 


7 




up to 25° 


9.0 


788 


1691 


88 


188 


12 


2 


25° - 40° 


3.7 


235 


718 


63 


194 


6 




40° - 50° 


1.5 


28 


324 


19 


216 


1 




up to 25° 


8.5 


694 


1584 


81 


186 


11 


3 


25° - 40° 


7.7 


372 


1533 


48 


199 


6 




40° - 50° 


6.3 


197 


1180 


31 


187 


4 



* - signal-to-noise ratio, S/N = N 7 (N 7 + 2 • N^ 1 / 2 , N 7 = ON - OFF, N cr = OFF 
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gle, for 7-ray showers has been inferred from Monte Carlo simulations as described in 
Konopelko et al., (1998). The energy threshold of the telescope system, defined as the 
energy at which the 7-ray detection rate reaches its maximum for the differential spec- 
trum dN 7 /dE ~ E -2 - 5 , is ~ 0.5 TeV at small zenith angles, and increases up to ~ 2 
TeV at 50° zenith angle. The rms error of the energy determination is AE/E ~ 0.18. 

In the stereoscopic observations the 
impact distance of the shower axis 
to a system telescope can be mea- 
sured with an accuracy < 10 m. 
The energy of a 7-ray shower is 
defined by interpolation over the 
"size" parameter S (total number 
of photoelectrons in image) at the 
fixed impact distance R, as E = 
fMc(S,R, 6), where 9 is the zenith 
angle and fMc is a function ob- 
tained from Monte Carlo simula- 
tions. Note that the Monte Carlo 
simulations include the sampling of 
detector response in detail. 
The energy distribution for the ON- 
and OFF-source events, after the 
orientation and shape image cuts, 
were histogramed over the energy 




E [TeV] 



Figure 1: The energy spectrum of Crab Nebula 
(preliminary) . 

range from 500 GeV to 30 TeV with 10 bins per decade. The 7-ray energy spectrum 
was obtained by subtracting ON- and OFF-histograms and dividing the resulting 
energy distribution by the corresponding collection area and the 7-ray acceptance. 
The Crab energy spectrum was then estimated separately for three observational pe- 
riods listed in Table 1. Note that Monte Carlo simulations have been adjusted for 
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each of the periods. As a systematic check of the spectrum evaluation technique, we 
compared the energy spectra measured from 4 telescopes (period 1) and 3 telescopes 
(period 2) data as well as with the spectra evaluated for different zenith angles. These 
spectra are in a good agreement, which shows the robustness of the technique. The 
preliminary Crab energy spectrum, produced by gathering data from three observa- 
tional periods, is shown in Figure 1. The collection area for 7-rays rises very fast 
in the energy range near the energy threshold of the telescope system, which is 500 
GeV. Even slight variations of the trigger threshold could lead to noticeable system- 
atic changes in the predicted spectrum behaviour in the energy range of ~ 0.5 — 1 
TeV. In addition, note that the statistical significance of detected 7-rays above 20 
TeV is not very high (~ 1.5a, N 7 =12). Thus in the present analysis we fitted the 
Crab data in the energy range 1-20 TeV by a simple power law, which yields 



The first set of errors in eqn. (1) is statistical and the second systematic is estimated 
using the Monte Carlo simulations which also take into account the current uncer- 
tainties in absolute calibration of the detector. The present Crab flux at 1 TeV, as 
measured with the HEGRA IACT system, is higher by a factor of 2 than estimated 
earlier for the HEGRA CT2 data using a cosmic ray calibration technique (Konopelko 
et al, 1996). A more detailed treatment of the telescope hardware, for example opti- 
cal smearing, increases the CT2 flux to a value (Petry et al, 1996) quite comparable 
to the present result. 

4. Discussion 

We derived the spectrum of the Crab Nebula from HEGRA data taken from 1997 
September to 1998 March. The data can be fitted by simple power law over the energy 
range 1-20 TeV. It does not exclude a possible slight steepening of the energy spectral 
usually predicted by inverse Compton modeling of the TeV 7-ray emission. Even a 
simple power law fit of the Crab data over the energy range 500 GeV to 25 TeV yields 
a differential spectrum index of ~ 2.45, which indicates a flatter spectrum. However, 
the noted change of the energy spectrum slope is within the current statistical and 
systematic errors and the data overall energy range are consistent with a simple power 
law fit in 1-20 TeV given by eqn. (1). Our data matches well at about 20 TeV the 
recent data published by CANGAROO group (Tanimori et al, 1998), which extends 
up to 70 TeV. Our data are consistent with a flat power law index of ~ 2.5 beyond 10 
TeV, which is likely to be for TeV 7-rays from 7r°s, generated by protons accelerated 
at the termination shock of the pulsar wind. 
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